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Abstract. We define nabla right fractional sum and diff'erence and obtain a nabla integration by 
parts formula. Some properties of nabla and delta fractional sums and differences are obtained to 
derive dual identities between the nabla and delta ones and other Q-dual identities to relate left and 
right ones. These dual identities give the impression that the definition of the delta right fractional 
difference presented by Abdeljawad and Baleanu in [T] and |12| are more appropriate than those 
introduced by other authors. Hence, the delta integration by parts formula formulated in |12| and 
the nabla one presented in this article are consistent. 

Keywords: right (left) delta and nabla fractional sums, right (left) delta and nabla Riomann 
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1 Introduction 

During the last two decades, due to its widespread applications in different fields of science and 
engineering, fractional calculus has attracted the attention of many researchers [9] 1101 1111 1231 1241 125| . 

Starting from the idea of discretizing the Cauchy integral formula. Miller and Ross [8] and Gray 
and Zhang ^18^ obtained discrete versions of left type fractional integrals and derivatives, called 
fractional sums and differences. Fifteen years later, several authors started to deal with discrete 
fractional calculus [1.3 4 5 ,6 f4l 1121 l20l 1211 [22l . benefiting from the theory of time scales originated 
by Hilger in 1988 (see il5j). 

In this article, we summarize some of the results mentioned in the above references and add 
more in the Caputo type and right type fractional differences. Throughout the article we almost 
agree with the previously presented definitions except for the definition of right fractional difference. 
We shall figure out that these definitions seem to be more convenient than the previously presented 
ones by proving some dual identities. These identities fall into two kinds. The first relate nabla type 
fractional differences and sums to delta ones. The second kind represented by the Q-operator relate 
left and right fractional sums and differences. Also, in the definition of left and right fractional sums 
we use the ones that agree with the general theory for time scales and not as in [5] or in |18| . This 
setting enables us to get identities resembling better the ordinary fractional case. Along with the 
previously mentioned points we are able to fit a reasonable nabla integration by parts formula which 
remains in accordance with the one obtained in 1121 but different from those obtained in [6] and 1141 . 

The article is organized as follows: The remaining part of this section contains summary to some 
of the basic notations and definitions in delta and nabla calculus. Section 2 contains the definitions 
in the frame of delta and nabla fractional sums and differences in the Riemann sense. Moreover some 
essential lemmas about the commutativity of the different fractional sum operators with the usual 
difference operators are established. These lemmas are vital to proceed in the next sections. The 
third section contains some dual identities relating nabla and delta fractional sums and differences 
in the left and right cases. Using these dual identities, power formulae for nabla left and right 
fractional sums and a commutative law for nabla left and right sums are obtained. In Section 4 an 
integration by parts formula for nabla fractional sums and differences is obtained with an accordance 
to the one obtained in 12^. Section 5 is devoted to delta and nabla Caputo fractional differences and 
their relations with the Riemann ones. Finally, Section 6 contains Caputo type fractional dynamical 
equations where a nonhomogoneous nabla Caputo fractional difference equation is solved to obtain 
nabla discrete versions for Mittag-Lefflcr functions. For the case a = I we obtain the discrete nabla 
exponential function 15^. In addition to this, the Q-operator is used to relate left and right fractional 
sums in the nabla and delta case. The Q-dual identities obtained in this section expose the validity 
of the definition of delta and nabla right fractional differences. 
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For a natural number n, the fractional polynomial is defined by, 

n-l 



t(") = r\it-j)= '■'^'> ^ , (1) 

where F denotes the special gamma function and the product is zero when t + 1 — j = for some j . 
More generally, for arbitrary a, define 

-r(t + i-a)' 

where the convention that division at pole yields zero. Given that the forward and backward differ- 
ence operators are defined by 

A/(i) = fit + 1) - fit) . SJfit) = fit) - fit - 1) (3) 

respectively, we define iteratively the operators A'" = A(A'"~^) and V" = V(V™~"'^), where m is 
a natural number. 

Here are some properties of the factorial function. 

Lemma 1.1. (Wi) Assume the following factorial functions are well defined. 

(i) At(") = 

(ii) it - At)t(f) = where G M. 
(Hi) /i^*^) = Vi^l + lY 

(iv) Ift<r, then i'"' < r'"' for any a > r. 

(v) IfO<a< 1, then > 

(vi) = it- f3)^°'h'^^). 

Also, for our purposes we list down the following two properties, the proofs of which are straight- 
forward. 

V4s-t)(°-i' =(a-l){p(s)-t)("-2). (4) 

Vt(p{.)-t)("-l) =-(a-l)(p(5)-f)(— 2). (5) 
For the sake of the nabla fractional calculus we have the following definition 

Definition 1.1. i^ fTgl flgl [T?! fTgf ) 

(i) For a natural number m, the m rising (ascending) factorial oft is defined by 

m— 1 

= Y[ it-k), t^ =1. (6) 
fc=0 

(ii) For any real number the a rising function is defined by 

r = I^i±i^, telR- {...,-2,-1,0}, 0"=0 (7) 

Regarding the rising factorial function we observe the following; 

(i) _ 

V(t") = at"-i (8) 

(ii) 

(P) = (t -Ha -!)(-). (9) 

(iii) _ 

At(.-p(i))" =-a(s-p{t))— 1 (10) 

Notation: 

(i) For a real o > 0, we set n = [a] + 1, where [a] is the greatest integer less than a. 
(ii) For real numbers a and b, we denote Na = {a, a -f 1, ...} and = {6, 6 — 1, ...}. 
(iii) For n S N and real a, we denote 

aA"/(i) ^ (-l)"A"/(t). 

iiv) For n G N and real b, we denote 

V?/{t) ^ (-l)"V"/(t). 
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2 Definitions and essential lemmas 

Definition 2.1. Let cr(t) = t + 1 and pit) = t — 1 be the forward and backward jumping operators, 
respectively. Then 

(i) The (delta) left fractional sum of order o > (starting from a) is defined by: 

_l t — a 

(a) The (delta) right fractional sum of order a > (ending at b) is defined by: 
b b 



(Hi) The (nabla) left fractional sum of order o > (starting from a) is defined by: 



v-"/W = -^ it-pi^)r-'fi^), teN^+i. (13) 

(iv)The (nabla) right fractional sum of order a > (ending at b) is defined by: 

6-1 6-1 

6V-"/(t) = — ^(s-p(t))-l/(«) = ^EW^)~*)°"'/W' *e (14) 

r(°) T^t r('^) 

Regarding the delta left fractional sum we observe the following; 

i) Aa " maps functions defined on Na to functions defined on Nn+a- 

ii) u{t) = Aq "/(t), n G N, satisfies the initial value problem 

A"u(t)=/(t), t e Na, u{a + j -1) =0, j = l,2,...,n. (15) 

iii) The Cauchy function — — vanishes at s = i — (n — 1), ...,t — 1. 

Regarding the delta right fractional sum we observe the following: 

i) i,A~" maps functions defined on to functions defined on 6-aI^- 

ii) u{t) = f,A~"/(i), n £ N, satisfies the initial value problem 

VXt) = /(t), te bN, uib-j + l) = 0, j = l,2,...,n. (16) 

iii) the Cauchy function ^''^"'j^j-i^^iy^ vanishes at s = t + 1, t + 2, t + (n — 1). 

Regarding the nabla left fractional sum we observe the following: 

i) Va " maps functions defined on N^, to functions defined on N^+i- 

ii) Va "/(t) satisfies the n-th order discrete initial value problem 

V"|/(t) = /(t), V'y{a)=0, i = 0,l,...,n-l (17) 

iii) The Cauchy function satisfies V"j/(t) = 0. 

Regarding the nabla right fractional sum we observe the following: 

i) i,V~" maps functions defined on [,N to functions defined on 6-lf^- 

ii) 6V~"/(t) satisfies the n-th order discrete initial value problem 

aA"y{t) = fit), aA'yia) = 0, i = 0, 1, n - 1. (18) 
The proof can be done inductively. Namely, assuming it is true for n, we have 

aA"+i 6V-("+i'/(i) = aA"[-A 6V-("+i'/(t)]- (19) 
By the help of lITOt . it follows that 

.A"+i 6V-("+i)/(t) = .A" 6V-"/(t) = fit). (20) 
The other part is clear by using the convention that X]fc=s = Oi s > t. 
(iii) The Cauchy function ^^~r\*y ^ satisfies aA"y(i) = 0. 
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Definition 2.2. (i)l^ The (delta) left fractional sum of order a > (starting from a ) is defined 
by: 

A^/W = A"A-("-"VW = r E teN,+(„_„) (21) 

rin — a) ^ — ' 

^ ^ s — a 

(ii) fist The (delta) right fractional difference of order q > (ending at b ) is defined by: 

1 

6A"/(i) = (-1)"V" 6A-("-")/(i) = J2 te 5_(„_,)N 

(22) 

( Hi ) The (nabla) left fractional difference of order a > ( starting from a ) is defined by: 



v°/(t) = v"v-<"-"Vw = ^7 r 5^ (t-p(s))"-"-VW, ten^+i (23) 

(iv) The (nabla) right fractional difference of order « > (ending at b ) is defined by: 

1 6-1 

bV^/W = (-1)"A" 6V-("-°'/(t) = ^ E(s - p(t))"-"- VW, t e i,_iN (24) 

Regarding the domains of the fractional type differences we observe; 

(i) The delta left fractional difference A^ maps functions defined on Na to functions defined on 

(ii) The delta right fractional difference (,A" maps functions defined on [,N to functions defined 
on ft_(„_„)N. 

(iii) The nabla left fractional difference V° maps functions defined on Na to functions defined 

on Na + l-n- 

(iv) The nabla right fractional difference (,V" maps functions defined on [,N to functions defined 
on b-i+„N. 

Lemma 2.1. J^' i^**!/ ct > 0, the following equality holds: 
A-"A/(i) = AA-"/(t) - 

r(a) 

Lemma 2.2. For any q > 0, the following equality holds: 

iA— Vb/(i) = Vt 6A-"/(i) - (^^^1^/(6). 

r(a) 

Lemma 2.3. JT^ For any « > 0, the following equality holds: 

V:^iV/(t) = - li^^±ll^/(a) (25) 

The result of Lemma l2.3l was obtained in [7] by applying the nabla left fractional sum starting 
from a not from a + 1. Next will provide the version of Lemma l2.3l bv applying the definition in this 
article. 

Lemma 2.4. For any o > 0, the following equality holds: 

V--V/(t) = VV-" - ^^'f^ ' f{a). (26) 
r(a) 

Proof. By the help of the following by parts identity 

V.[(t - s)~fis)] = V4i - s)~f{s) + {t- p(s))~V./(s) (27) 
= ^(a - l){t - p{s)r^fis) + (t- p{s)Y^^sf is) 



we have 



v-"v/(t) = -^ J2 it-p{s)r-'^sf{s) 

r(o) ^ 



s — a-j-1 



+ («_!) j2 {t-Pis)r-^m]= (28) 

r(a) 



s — a + 1 



(t - a)"-! 1 



^ ' ^ ' s— a+1 

On the other hand 



vv--/{t) = 

t;^ E Vt(t-p(s))~/(s)= \ Yl (30) 

□ 

Remark 2.1. Let a > and n = [a] + 1. Then, by the help of Lemma \2.4\ we have 

^^afit) = VV"(V-("-"'/(t)) = V"(VV-("-"'/(t)). (31) 

or 

VV^/(i) = V"[V^("-"'v/(t) + ^— ^ —f(a)] (32) 

r(n — a) 



Then, using the identity 



(t-a)""--^ _ (t-a) — 1 



r(n-a) r(-o) 

lue m/er 4/ia4 f26|) is valid for any real a . 



(33) 



By the help of Lemma [Zil Remark |2?T] and the identity V(i - a)""^ = {a - l){t - a)"''^ , vie 
arrive inductively at the following generalization. 

Theorem 2.5. For any real number a and any positive integer p, the following equality holds: 

V-'^ ^/(i) = V''V--/(t) - J2 FT^^T— ?TT^'-^('^)- (34) 

^r(a + fc-p+l) 

where f is defined on Na ■ 

Lemma 2.6. For any o > 0, the following equality holds: 

aAfit) = aA 6V-"/(t) - i^^l^/(fe) (35) 

r(a) 

Proof. By the help of the following discrete by parts formula: 

AsMs)-p{t))~f{s)] = 
(a - l){s - p[t)r^f{s) + {s- p{t)r^Af{s) (36) 

1 6-1 

iV-" aAfit) = -— E(. - p(t))"-i A/(.) = 

b-l b-l 

- E A,((p(s) - p{t)r~^f{s)) + (a - 1) E(. - p(t))"-2/(s)] = (37) 



we have 



r(")' 



s = t 



b-l 



— L— y(s-p(t))— 2/(,)-i^^ll-^/(fe). (38) 



On the other hand, 

.A i,V-"/(t) 



b-l 6-1 

r(a) ^ r(a - 1) ^ 

where the identity 

At (s - p(t))~ = -(a - l)(s - p(t))~ 
and the convention that (0)"~^ = are used. □ 
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Remark 2.2. Let a > and n = [a] + 1. Then, by the help of Lemma \2. 6\ we can have 

aA b^'^fit) = aA aA"i bV-'^-^'/W) = aA"( aA 6 V" /(* )) (40) 

or 

(h — V-^Ti — a — 1 

„A i,V"/(t) = aA"[ i,V-("-°' aAfit) + —f{b)] (41) 

r(n — a) 

Then, using the identity 



„A"-i = -i^ (42) 

r(n-a) r(-a) ^ ' 

we infer that 1 35(1 is valid for any real a. 

By the help of Lemma [T6l Remark [231 and the identity A(fe - t)""! = -{a - 1)(6 - i)°~^, if 
wc follow inductively we arrive at the following generalization 

Theorem 2.7. For any real number a and any positive integer p, the following equality holds: 

P^l /I _ ,\a-p + k 

i,V- aAPfit) = „Af i,V-"/(i) - J2 FT^T VTT "^'/W (^3) 

^r(a + A:-p+l) 

where f is defined on ijN and we remind that aA'=f{t) = {-!)'' A'' f{t). 

3 Dual identities for right fractional sums and dif- 
ferences 

The dual relations for left fractional sums and differences were investigated in Indeed, the 

following two lemmas are dual relations between the delta left fractional sums (differences) and the 
nabla left fractional sums (differences). 

Lemma 3.1. Let < n~l < a < n and let y{t) be defined on Na- Then the following statements 
are valid. 

(iKA<S)y{t -a) = V^yit) for t e 

(ii) {A-°')y{t + q) = V^"j/(t) for t G N,. 

Lemma 3.2. {Sf Let 0<n — f<Q!<n and let y(t) be defined on 'Ha—n- Then the following 
statements are valid. 

(i) Al_.^y(t) = (V«_„j/)(t + a) for t G N_„. 

(ii) A-JX"\{t) = (V-i"„-"'j/)(« -n + a) fort G No. 

We remind that the above two dual lemmas for left fractional sums and differences were obtained 
when the nabla left fractional sum was defined by 

^a^f(t) = -^i2(^-p{s))~f{s), t€Na (44) 

Now, in analogous to Lemma 13. f I and Lemma 13.21 for the right fractional summations and 
differences we obtain 

Lemma 3.3. Let y{t) be defined on j,+iN. Then the following statements are valid. 

(i) itA")y{t + a)= t+i^"y{t) forte 6_„N. 

(ii) { tA-°')y{t - a) = b+iV-"s/(i) /ort G ^N. 

Proof. We prove only (i). The proof of (ii) is similar and easier. 

( i,A")y(t + a) = i-irw" 6 A-("-")y(t + a) = 

^ ^ (,_j„l_c.)("-"-l)y(s) = ^-i ^^(s-t-l-+n-a)("— l)y(.) (45) 

^ ' s — t-\-n 
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Using the identity = (t -\- a — , we arrive at 



( ,A'^)y{t + a) = ^(s - p(t))"— = 

r(n -a) fri 

(-1)"A" 6+iV-("-°'2/(t) = 5+iV"j/W. (46) 

□ 

Lemma 3.4. Let 0<n — l<a<n and let y{t) be defined on Then the following 

statements are valid, 
(i) 

(ii) 

„_c.A-("-"'j/(i) = „_^+iV-("-°)y(t + n-«), te oN 
Proof. We prove (i), the proof of (ii) is similar. By the definition of right nabla difference we have 



b 



n-a. + l^"y{t-a)= „A"- J2 (s-p{t-"))"-"-'yW = 

^ n — CL _| n — (H 

aA"- Y. (s-p(t-")r-"-^S/W = V?- ^ (,_p(t + n_„))"-"-iy(s) 

Fin — ol) ^ T{n — a) ^ 

^ ' s—t — a ^ ' s — t-\-n — a 

(47) 

By using ^ it follows that 

^ n — a 

„_,+iV"j/(t-a) = V?- J2 (s-'T{t)Y"-''-'^y{s)= n-cA^y{t) (48) 

r(n - a) f--' 

□ 

Note that the above two dual lemmas for right fractional differences can not be obtained if we 
apply the definition of the delta right fractional difference introduced in If 41 and |6] . 

Lemma 3.5. HH' Let a > 0, ^ > 0. Then, 

i,_^A-"(fe - t)(^) = ^^!i±2L^{b - i)(^+") (49) 
r(fj. + a + l) 

The following commutative property for delta right fractional sums is Theorem 9 in If2| . 
Theorem 3.6. Let o > 0, fi > 0. Then, for all t such that t = b — {fi + a) [mod 1), we have 

i,A-"[ 6A-''/(t)] = bA-(^+")/(t) = i,A-''[bA--/(t)] (50) 
where f is defined on f^N . 

Proposition 3.7. Let f be a real valued function defined on 5N, and let a,(5 > 0. Then 

bV-"[ 6V-'3/(t)] = ,V-("+«/(t)= 6V-'9[,V--/(i)] (5f) 
Proof. The proof follows by applying Lemma l3.3f ii) and Theorem 13.61 above, fndeed, 
i,V— [ bV-'9/(t)] = b-i^-^f(t -13) = 

6-1 A-" 6-1 A-'3/(t - (a + 13)) = A-(-+«/(t - (a + ^t)) = t V-(°+«j/(t) (52) 

□ 

The following power rule for nabla right fractional differences plays an important rule. 
Proposition 3.8. Let a, fi > 0. Then, for t £ , «ie have 

,V-"(fe - t)-p = -I(^±il-(6 - t)^ (53) 
r(a + fj. + l) 
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Proof. By the dual formula (ii) of Lemma 13. 31 we have 

b^-°'{b-tf= t-l^-"{b-rf\r^t-c = 
1 ''-^ 

Then by the identity = (t + a — and using the change of variable r = s — + 1, it follows 

that 

tV-'-(b-tf = 

b — fi 



-L_ J2 (r-a(t-a-/. + l))(--i)(6-r)" = (b_^A-"(6-«)'')U=_,_^+i+t. (55) 
Which by Lemma 13. 51 leads to 



r(M + i) 



W ^ (fe-t + « + ;.-l)("+^) = ^(^ + ^) (fc-f)"+M (56) 

r(a + /i + i) r(Q + /j + i) 

□ 

Similarly, for the nabla left fractional sum we can have the following power formula and exponent 

law 

Proposition 3.9. Let a, fj, > 0. Then, for t £ Na , we have 

(t - af = {t - a)~ (57) 

r(o + ^ + 1) 

Proposition 3.10. Let f be a real valued function defined on Na, and let a, [5 > 0. Then 

V-"[V-'5/(t)] = V-("+«/(t) = V-'3[V-"/(t)] (58) 

Proof. The proof can be achieved as in Theorem 2.1 [5], by expressing the left hand side of JSSj, 
interchanging the order of summation and using the power formula l|57| l. Alternatively, the proof 
can be done by following as in the proof of Proposition 13. 71 with the help of the dual formula for left 
fractional sum in Lemma 13. II after its arrangement according to our definitions. □ 

4 Integration by parts for fractional sums and dif- 
ferences 

We first state the integration by parts for delta fractional sums and differences. 

Proposition 4.1. f jg/ Let a > 0, a, b € R such that a < b and b = a + a {mod 1). // / is defined 
on Na and g is defined on i,N, then we have 

b b—a 

J2 (A-«/)(s)g(.) = ^/(s)6A-"9(.). (59) 



s — a + a 



Proposition 4.2. flSf Let a > be non-integer and assume that b = a + {n — a) (mod 1). // / is 
defined on bN and g is defined on Na, then 

b-n + l b-(?i-a) + l 

J2 f{s)A^g{s)= Yl 9is) b^'^fis). (60) 

s — a + {n — a) — 1 s — a + n— 1 

Proposition 4.3. For a > 0, a, b 6 R, / defined on Na and g defined on jN, we have 

6-1 6-1 

E 9WV^"/W= E /{^)6V'-9(s). (61) 
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Proof. By the definition of the nabla left fractional sum we have 

b-l -. 6-1 s 

E f = E 9{s) J2 («-pwr"VM. (62) 

s — a + l ^ ^ s — a + 1 r — a + 1 

If we interchange the order of summation we reach at ( I61I I . 

□ 

By the help of Theorem [231 Proposition ISTTOl l ITTll and that Va '"""'/(a) = 0, we can, for the 
nabla left sums and differences, obtain 

Proposition 4.4. For a > 0, and f defined in a suitable domain Na, we have 

^a^a'^fit) = fit), (63) 

V-"V?/(t) = /(t), whena<^N, (64) 



and 

n-l , 

V-°'\7"f(t) = f(t) - V ^ 

k\ 



V-°V^/(t) = fit) - J2 ^^Tr-VV(a), .when a = n g N. (65) 



We recall that II64I I is valid in the usual fractional case for sufficiently good functions such as 
continuous functions. As a result of this it was possible to obtain an integration by parts in the 
Riemann fractional derivative case for certain class of functions (see page 76, and for more 
details see [TO]). Since discrete functions are continuous we see that the term f{t)\t=a, 
for < o < 1 disappears in I I64II . with the application of the convention that X]s=a+l /(^) ~ ^■ 
However, in the Caputo case the initial type conditions starting from a appear as we shall see in the 
next sections. By this connection, we would like to ask the reader to compare with 0, where the 

term Va "V(*)|t=£i = f{<^) appears and hence it was possible to employ the Riemann type initial 
value fractional difference equation to obtain a discrete fractional version of Gronwall's inequality. 
However, we remind that the authors there obey | 18] in defining the nabla left fractional sum. In our 
case the Caputo fractional difference will be the more suitable tool to obtain such fractional version 
of Gronwall's inequality. 

By the help of Theorem [TtI Proposition |321 iO and that t,V-("-"'/(b) = 0, we can, for the 
nabla right sums and differences, obtain 

Proposition 4.5. For a > 0, and f defined in a suitable domain jN, we have 

i,V-"/(t) =/(t), (66) 

bV""" 6V"/(t) = fit), when Q ^ N, (67) 

and 

(I, _ f)k 

bV-°' 6V"/(t) = fit) - J2 a^'' fib), when a = n e N. (68) 

fc=0 

Proposition 4.6. Let a > be non-integer. If f is defined on i,N and g is defined on Na , then 

6-1 6-1 

E /WVaS(«)= E 9is)t^"fis). (69) 

s — a+l s — a+l 

Proof. By the help of equation I I67I I of Proposition 14. 51 we can write 

6-1 6-1 

E /(s)A"g(s)= E 6V-( 6V"/(s))V^g(s) (70) 

s — a + l s— a + l 

and by Proposition 14.11 we have 

6-1 6-1 

E /(s)A-g(s)= E 6V"/(s)V--V-s(s). (71) 
Then the result follows by equation (|64| l of Proposition 14.41 □ 
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5 Caputo fractional differences 



In analogous to the usual fractional calculus we can formulate the following definition 

Definition 5.1. Let a > 0, a ^ N. Then, 

the delta a—order Caputo left fractional difference of a function f defined on Na is defined 

by 

^A^/(t) ^ A-("-"'A"/(t) = (t-a{s))(^-'^-^^A^f{s) (72) 

r(n - a) 

(ii) [Ijj the delta a— order Caputo right fractional difference of a function f defined on (,N is 
defined by 

1 '' 

fA-ZW^ ,A-("— )V?/(i) = J2 («-'^(i))<""°"'^V?/(s) (73) 

Fin — a) ^—^ 

where n = [o] + 1. 

// a = n e N, then 

^A-/(t) ^ A"/(t) and f A"/(t) 4 V^/{t) 

(iii) the nabla o— order Caputo left fractional difference of a function / defined on Na is defined 

by 

t-(n-a) 

(iv) the nabla o— order Caputo left fractional difference of a function / defined on {,N is defined 

by 

1 b-l 

f V-/(t) ^ ,V-("— ) aA"/(i) = ^(s - p(i))"-"-i aA"f{s) (75) 

r(n - ") jr^ 

If a = n e N, then 

^V-/(i) 4 V"/(i) and f V"/(t) ^ aA"/(t) 

It is clear that "-^A" maps functions defined on Na to functions defined on N;,^(„_q), and that 
pA" maps functions defined on to functions defined on [,_(„_q)N. Also, it is clear that the nabla 
left fractional difference V" maps functions defined on Na to functions defined on Na+i_„ and the 
nabla right fractional difference i,V°' maps functions defined on to functions defined on 5_i+„N. 

Riemann and Caputo delta fractional differences are related by the following theorem 



Theorem 5.1. JYj' For any a > 0, we have 



n-l 



CA-/(i) = A^/(t) - J2 r(fc "r+i) ^'-^^"^ ^^^^ 



and 



"(fc-a + 1) 



i-i (I _ i-)(,k-a) 

^A'^fit) = ,A"/(t) - J2 ^u ^ifC^y (77) 

fe=0 a + i) 

In particular, when < o < 1, we have 



Aafit) = A2f{t) - ' ' f(a). (78) 
r(l - a) 

^Af{t) = ,A-/(i) - ^±zJ}lZ^f(b) (79) 
r(i - a) 

One can note that the Riemann and Caputo fractional differences, for < « < 1, coincide when 
/ vanishes at the end points. 

The following identity is useful to transform delta type Caputo fractional difference equations 
into fractional summations. 
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Proposition 5.2. ^ Assume a > and f is defined on suitable domains Na and ;,N. Then 

-OL CaO. 



A 
and 

6_(„_„)A- fA"/(t) = fit) - J2 J Vg/(fe) (81) 

/ra particular, if < a < 1 then 

^;+(n-.) ''^a/W = fit) - fi^) '^nd 6_(„_,) A— f A-/(t) = /(t) - fib). (82) 
Similar to what we have above, for the nabla fractional differences we obtain 
Theorem 5.3. For any a > 0, we have 

^V-/(i) = V^/(t) - J2 p _,-,, VV(a) (83) 

^^r(fc-a + l) 



and 



'V^fit) = ,V"/(i) - ^ aAVC'). (84) 

^^Tik-a + 1) 



In particular, when < q < 1, we have 



and 



"^afit) = V^fit) - " fia) (85) 

r(i - a) 

C'V/(i)= ftV"/(f) - ; ° /(b) (86) 

r(i - o) 

Proof. The proof follows by replacing a by n — a and p by n in Theorem 12.51 and Theorem 12.71 
respectively. □ 

One can see that the nabla Riemann and Caputo fractional differences, for < a < 1, coincide 
when / vanishes at the end points. 

Proposition 5.4. Assume q > and f is defined on suitable domains Na and [,N. Then 



it -a) 

^ a aJ y-j = J y-j - / , 

k = 

and 



^V^/(i) = fit) - ^ ^^vV(a) (87) 



f V-/(t) = fit) - ^ ^iL , A*/(b). (88) 

In particular, if < a < I then 

'^V^/(t) = fit) - fia) and .V"" f V"/(t) = /(i) - /(fe) (89) 

Proof. The proof of II87I I follows by the definition and applying Proposition 13.101 and II65I I of Propo- 
sition 14.41 The proof of l|88|l follows by the definition and applying Proposition 13.71 and H68|l of 
Proposition 14.51 □ 

Using the definition and Proposition [3]9] and ProDosition l3.8l we can find the nabla type Caputo 
fractional differences for certain power functions. For example, for 1 7^ /3 > and a > we have 

C (i - a)^ = -^^^ it - a)l^^^ (90) 

and 

C V- (fe - t)'^ = (b _ t)/^^^. (91) 



However, 
whereas 



fV°l = (92) 

r(l-a) r(l-a) ^ ^ 

In the above formulae 1 190 I I and II91I I. we apply the convention that dividing over a pole leads to 
zero. Therefore the fractional difference when /3 — 1 = a — j, j = 1, 2, n is zero. 
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Remark 5.1. The results obtained in Theorem 15. Ib nd afterward agree with those in the usual 
continuous case (See f77| pages 91,96). 



6 The Q-operator and fractional difference equa- 
tions 

If f{s) is defined on Na n i,Af and a = b {mod 1) then (Q/)(s) = f{a + b — s). The Q-operator 
generates a dual identity by which the left type and the right type fractional sums and differences 
are related. Using the change of variable Ji = a-|-b — s, inlTlit was shown that 

A-"Q/(t) = Q ftA — /(t), (94) 

and hence 

^A?Q/(t) = Q(C^A"/)(t). (95) 
The proof of I I95I I follows by the definition, I I94I I and by noting that 

-QV/(i) = AQfit). 
Similarly, in the nabla case we have 

V-"Q/(t) = Q 6V-"/(t), (96) 

and hence 

c^V^Q/(t) = Q(fV"/)(i). (97) 
The proof of II97I I follows by the definition, II96I I and that 

-QAf{t) = \7Qf(t). 

The Q-dual identities I I95I I andf I97l l are still valid for the delta and nabla (Riemann) fractional 
differences, respectively. The proof is similar to the Caputo case above. 

It is worthwhile to mention that the Q-dual identity II95I I can not be obtained if the definition of 
the delta right fractional difference introduced by Nuno R.O. Bastos et al. in 1141 or by Atici F. et al. 
in [6]. Thus, the definition introduced in [T] and 1121 is more convenience . Analogously, the Q-dual 
identity I I97I I indicates that the nabla right Riemann and Caputo fractional differences presented in 
this article are also more convenient. 

It is clear from the above argument that, the Q-operator agrees with its continuous counterpart 
when applied to left and right fractional Riemann Integrals and the Caputo and Riemann deriva- 
tives. More generally, this discrete version of the Q-operator can be used to transform the discrete 
delay-type fractional functional difference dynamic equations to advanced ones. For details in the 
continuous counterparts see [2]. 

Example 6.1. [Tf Let 0<a<l, a = a — 1 and consider the left Caputo nonhomogeneous fractional 
difference equation 

'^Kvit) = Mj{t + a - 1) + fit), y{a) = ao, t G Nq. (98) 

Note that the solution y{t), if exists, is defined on Na and hence ^ A"y(t) becomes defined on 
T^a+(l—a) = Nq. Thus, if We apply Aq" on the equation i98V and then use the successive approxi- 
mation the following delta type solution is obtained 

yit) = aoE^^^{X,t) + Y, E(^^^^(X,t - a{s))f{s) (99) 

s=0 

where the delta type discrete Mittag-Lefffer functions are defined by 

Definition 6.2. (Delta Discrete Mittag-Leffter) 13' For A G K and a,l3,z eC with Re(a) > 0, the 
Delta discrete (like) Mittag- Leffier functions are defined by 

^0 T{a.k + I3) ^ ' 

For p = 1, it is written that 



12 



Next, we solve nabla Caputo type nonhomogeneous fractional difference equation to formulate 
nabla type discrete Mittag-LcfHer functions . Those functions generalize the nabla discrete exponen- 
tial functions. For details about delta and nabla type exponential functions on arbitrary time scales 
we refer the reader to f|15|, pages 10, 76). 

Definition 6.3. (Nabla Discrete Mittag-Leffler) For A € K and a,/3,z € C with Re{a) > 0, the 
nabla discrete (like ) Mittag-Leffler functions are defined by 



E-^{X, z) = J" X'' — -. (102) 

For j3 = 1, it is written that 

i?^(A,.)^i?^{A,.) = X:A'=^^^^. (103) 

Example 6.4. Let 0<a<l, aSM and consider the nabla left Caputo nonhomogeneous fractional 
difference equation 

^aVit) = ^y{t) + f(t), y(a) = ao, te Na. (104) 
// we apply Va " on equation fTOj}) then by iSM we see that 



y{t) = ao + -^ (t_p(s))— iy(s) + V-"/(t). 

To obtain an explicit solution, we apply the method of successive approximation. Set yo{t) = ao and 

ym(t) =ao + ^-"[Xym-lit) + /(*)], m = 1, 2, 3, .... 
For m = I, we have by the power formula {57 



m(t)=«o[l+ ^^* °ij + V-°/(t). 
r(a + 1) 

For m = 2, we also see by the help of Provosition \3.10\ that 

y2{t) = ao + AV-"[ao + ^*~°\" ] + V--/(t) + AV-2-/(t) = 
r{a + 1) 

Proceeding inductively and making use of Provosition VS . 10\ and let m oo we obtain the solution 

X'^(t - a)*" 
V{ka + 1) 



Then, 

t 



y{t)=aoE-{X,t^a) + f2^' , , X] (105) 

Interchanging the order of sums in \1U5\) . we conclude that 



yit) = aoE^^^iX,t-a)+ ^ ^ A" rfl^ \ (^O^) 



That is 

t 

y(t) = aoE^{X, t-a)+ ^ ^^^^(A, t - p{s))f{s). (107) 

s — a + l 

Remarlc 6.1. // we solve the nabla discrete fractional system {104'i with a = 1 and ao = 1 we 
obtain the solution 

^w-£a^^+ t EA^^^^/c.). 

fe = ■ s=a+lfc=0 

The first part of the solution is the nabla discrete exponential function e\{t, a). For the sake of 
more comparisons see f \15t . chapter 3). 
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Example 6.5. '1] Let 0<cj<l, a = a — I and b such that a = b (mod 1). Let y(t) be defined on 
Na n fcN. Consider the following Caputo right fractional difference equation 

bA"Qy{t) = \y{2a + b^t),+fia + b-t) {Qy)(b) = ao. (108) 

If we apply the Q operator on the Caputo fractinal difference equation n08\} . then we obtain the left 
Caputo fractional difference equation iQSV . For more about the Q-operator and its use in functional 
fractional differential equations we refer to ]2& . 

Similar to EaxamDle l6.5l above. we can use the Q-operator to transform, as well, nabia left type 
fractional difference equations to right ones and vise versa. 
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